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Abstract 

The peeling behaviour of the Weyl tensor near null infinity is determined for an asymptotically 
flat higher dimensional spacetime. The result is qualitatively different from the peeling property 
in 4d. To leading order, the Weyl tensor is type N. The first subleading term is type II. The 
j—i next term is algebraically general in 6 or more dimensions but in 5 dimensions another type N 

I ■ term appears before the algebraically general term. The Bondi energy flux is written in terms of 

"Newman-Penrose" Weyl components. 

^ ! 1 Introduction 

CO. 

t~^. . In 4 spacetime dimensions, the Weyl tensor of an asymptotically flat spacetime exhibits the "peeling" 
property: 

<W = X-'C^ + A- 2 cf// CT ) + \-*CjL% + X-^X + 0(A" 5 ) (1) 



where indices fj,,v,... refer to a basis parallely transported along an outgoing null geodesic with affine 
parameter A. In the first term, C^\ a is a Weyl tensor of algebraic type N and the subsequent terms 
involve Weyl tensors of algebraic types III, II and I. The tangent to the geodesies is the repeated 
principal null direction for the type N, III and II terms, and a principal null direction for the type I 
term. 

This result was originally derived using Bondi coordinates [U[2]- In this approach, one assumes 
that the metric components can be expanded in inverse powers of a coordinate r. A more geometrical 
proof can be given using the definition of asymptotic flatness in terms of a conformal compactification 
OH]. In this case, the result follows from the assumed smoothness of the unphysical spacetime. 
This smoothness assumption (or the assumption of an expansion in inverse powers of r) excludes 
some spacetimes in which radiation is present near spatial infinity. In this case, the peeling property 
is modified by an 0(A _4 logA) term [HE] (see also |7J)Q 

In d > 4 dimensions, a definition of asymptotic flatness at null infinity using conformal com- 
pactification is possible only for even ci [El |9]. It has been shown that this definition is preserved 
by linearized metric perturbations arising from compactly supported initial data [8]. It has also 
been argued that a vacuum spacetime satisfying this definition arises from initial data describing 
a small (but finite) perturbation of Minkowski spacetime that coincides with Schwarzschild initial 
data outside some compact set |10| . Just as in 4d, there are more general initial data that do not 
give a smooth null infinity |llj . 



1 We thank M. Dafermos for pointing out these references. 
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For odd d, conformal compactification is unsatisfactory because the unphysical spacetime cannot 
be smooth in any radiating spacetime [12]. Instead, one can follow the older approach of defining 
asymptotic flatness at null infinity using Bondi coordinates [131 E] . In section [2] we will weaken this 
definition slightly and demonstrate equivalence of the conformal and Bondi definitions for even d. 
For odd d, it remains to be shown that there exists an interesting class of initial data that gives rise 
to a spacetime which satisfies this definition. 

The goal of this paper is to determine how the Weyl tensor peels near null infinity in a spacetime 
satisfying one of the above definitions of asymptotic flatness. As just mentioned, at least for even 
d, this includes a large class of physically interesting spacetimes, but probably also excludes some 
physically interesting spacetimes. However, we can hope that in the latter case, just as in 4d, the 
peeling behaviour is modified only at a sufficiently high order that our result is still useful. 

Two previous papers have investigated peeling using the conformal approach to asymptotic 
flatness [15\ I16j. Both papers concluded that peeling is similar to the d = 4 case. They started from 
the assumption that all components of the unphysical Weyl tensor decay at the same (unspecified) 
rate near null infinity. However, Ref. [12J showed that this assumption is not true even for linearized 
perturbations of Minkowski spacetime, and argued that peeling should be qualitatively different for 
d > 4. This is what we find. 

In section [3j we determine the behaviour of the Weyl tensor near null infinity in a spacetime 
satisfying the "Bondi definition" of asymptotic flatness (since this is valid for odd or even d and 
equivalent to the conformal definition for even d). For d > 6, we find the following result: 

<W = A-t^cM, + \~ d/2 C«% + X-^C^l + .... (2) 

Again A is an affine parameter along a null geodesic and refer to a parallelly transported 

basis. The superscripts N,II,G refer to the higher-dimensional classification of the Weyl tensor 
developed in Ref. [IT] , based on the concept of Weyl Aligned Null Directions (WANDs). A type N 
or type II Weyl tensor admits a "multiple WAND", in this case it is the tangent to the geodesic. 
The type II term in ([2]) is not the most general type II Weyl tensor: it obeys additional restrictions 
explained below. Type G denotes an algebraically general Weyl tensor. The ellipsis in ([2]) denotes 
terms of order \-W 2 + 2 ) (even d) or \-W 2 +V 2 ) (odd d). 

For even d, the derivation of this result requires no more than the definition of asymptotic 
flatness. We do not use the Einstein equation, so this result is valid for any energy-momentum 
tensor consistent with asymptotic flatness. For odd d, we need to use some additional information 
from the Einstein equation: a mild condition on the decay of the Ricci tensor near null infinity is 
required to eliminate a term of order \~( d l 2+l l 2 ) from ([2]). 

The case d = 5 is exceptional. In this case, the Einstein equation no longer eliminates the term 
of order \~( d l 2+1 / 2 ) = A~ 3 . Instead, it fixes this term to be quadratic in the leading order metric 
perturbation and hence non-zero in any radiating spacetime. The result is that an additional type 
N term appears between the type II and type G terms: 

CW = A" 3 / 2 * + + A- 3 CS; + X-y^U + 0{\-% (3) 

The subleading type N term is distinct from the leading order type N term. The presence of this 
term can be attributed to the nonlinearity of the Einstein equation. For d > 5, nonlinear effects 
decay faster and this term does not arise. 

Refs. [U [9] [13] [T4] gave expressions for the rate of change of the Bondi energy at null infinity. 
For applications (e.g. higher-dimensional numerical relativity) it is convenient to have results that 
can be calculated easily and do not refer to a particular coordinate chart. This can be achieved by 
writing the result in terms of the asymptotic Weyl tensor components. We do this in section 0] 
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2 Definitions of asymptotic flatness 



2.1 Conformal definition for even d 

For even d > 4, Ref. [8] defined a spacetime (M,g) to be asymptotically flat at null infinity as 
follows. Given the (physical) metric g and the Minkowski metric rj, we would like to specify the 
precise rate at which g approaches 77 asymptotically. We do this by conformally compactifying both 
M and Minkowski spacetime so that "infinity" is now at a finite metric distance. Thus, we obtain 
the "unphysical" spacetime (M,g) and the "background" spacetime (M,g), where the metrics g and 
g are related to the respective physical and flat metrics via 

g ab = 9? g ao , 9ab = ^ 2 r] a b (4) 

with the conformal factor Q 2 satisfying the usual suitable properties. 

Now, the spacetime is defined to be asymptotically flat at null infinity if 

9ab - 9ab = O^ 2 - 1 ), l ai ... ad - e ai „^ = 0(^ 2 ), 

(g ab - g ab )m) a = 0(^/ 2 ), (g ab ~ g ab )(dn) a (dn) b = 0{^l d l 2+1 ), (5) 

where g ab and g ab are the inverse metrics of g and g, respectively and e and e are the volume forms 
on (M,g) and (M,g), respectively. Following Ref. [4], if L a b...c is a tensor field on M then the 
notation -L a 6... c = (D(£l s ) means that Q~ s L a b^ c is smooth at future null infinity. 

2.2 Definition using Bondi coordinates 

For general d > 4, Ref. |14| defined a spacetime to be asymptotically flat at future null infinity if, 
outside some cylindrical world tube, coordinates (u, r, x 1 ) can be introduced following [2 J such that 
the metric takes the form 

ds 2 = -Ae B du 2 - 2e B dudr + r 2 /i /J (dx / + C 1 du)(dx J + C J du) (6) 

with 

det hjj = det w/ j (7) 

where ujj(x) is the unit round metric on S d ~ 2 . Surfaces of constant u are null with topology 
K x S d ~ 2 where x 1 are coordinates on S d ~ 2 and R corresponds to the null geodesies generators of 
the surface. These generators are given by u, x 1 = constant and r is a (non-affine) parameter along 
the generators. A, B, C 1 and hjj are functions of all of the coordinates. It is assumed that, at large 
r, they be expanded in inverse powers of r (even d) or y/r (odd d) withH 

A = l + 0(r-W 2 -V), B = 0{r- d ' 2 ), 
C 1 = 0(r- d ^ 2 ), h I j=cj I j + 0(r-W 2 - 1 ^). (8) 

For odd d it appears that an extra condition is required (discussed for d = 5 in Ref. p3]). One way 
of seeing this is to note that the results of Refs. [U [12] suggest that, for linearized perturbations of 
Minkowski spacetime (arising from compactly supported initial data), the components of the metric 
perturbation each will be some half- integer power of 1/r times a smooth function of 1/r. Hence 
each component will involve either integer powers of 1/r or half-odd-integer powers, but not both. 
Therefore, the presence of both integer and half-odd-integer powers in the expansions of individual 
metric components can be attributed to nonlinear effects. One would expect these only to affect 
terms beyond a certain order in the above expansions. If so, at low enough order, these expansions 
should contain only integer powers, or only half-odd-integer powers. This is indeed the case if one 
imposes the additional boundary condition that the expansion of hu in inverse powers of \[r contains 
no term of order r - ^/ 2-1 / 2 ) (see below). 

2 Ref. [H] took B = 0(r- d ), which was obtained by solving the vacuum Einstein equation. We have weakened 
this condition since we don't want to assume vacuum. 
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2.3 Equivalence of definitions for even d 

Starting from the Bondi definition, define Q = 1/r. It is straightforward to show that this satisfies 
the conformal definition with (conformally flat) background metric 

g = -tftdu 2 + 2dudft + ui J dx 1 dx J . (9) 

Now consider a spacetime that is asymptotically flat according to the conformal definition. Write 
the flat metric in the form 

rj = -dU 2 - 2dUdR + R 2 uj IJ (X)dX I dX J . (10) 

Now define Q = 1/R. The background spacetime is 

g = n 2 7] = -n 2 dU 2 + 2dUdQ + uj IJ (X)dX I dX J . (11) 

X + is at $7 = and O > corresponds to the spacetime interior. 

In these coordinates, the definition of asymptotic flatness reduces to the following conditions on 
the unphysical spacetime 

guv = -n 2 + o(n d / 2+1 ), g un = i + o(n d / 2 ), g UI = o(n d / 2 ) 

gnn = 0(n d / 2 - 1 ), ~g m = O^ 2 " 1 ), ~g I j = u I j + 0(n d / 2 - 1 ) (12) 

and 

det gu = detuu + 0(fl d/2 ). (13) 

Now convert to Gaussian null coordinates based on the null surface X + in the unphysical spacetime 
as follows. Consider the (past-directed) null geodesic (of g) that passes through the point on X + 
with coordinates [u, 0, x 1 ) and has tangent vector d/dQ there. Let A denote the affine parameter 
along the geodesic. Since g is required to be smooth near I + it follows that the coordinates along 
the geodesic are smooth functions of A in a neighbourhood of A = 0. Expanding them in a Taylor 
series in A and substituting into the geodesic equations gives 

u = u + o(\ d/2 ), n = x + o(\ d/2+1 ), x 1 = X I + o(\ d/2 ). (u) 

We take (u, A,x 7 ) as new coordinates. In these coordinates, the unphysical metric is 

-A 2 + e>(A d/2+1 )] du 2 + 2dudX + OiX^dudx 1 
uu(x) + OiX^ 2 - 1 )] dx I dx J (15) 
where all components are smooth at A = and 

det g u = detu I j + 0(X d/2 ). (16) 
We now replace A with a non-affine parameter r defined by 

det gu N 1/(2(d ~ 2)) 

r = \l 1 

so 



+ 



^-i^^^j =A- 1 (l + 0(A d / 2 )j (17) 

X = r- 1 (l + 0(r- d / 2 )) (18) 



and 



1 =r (l + 0{r- d ' 2 )^ . (19) 

In coordinates (u, r, x 1 ) the physical metric takes the Bondi form (|6l7p with metric coefficients that 
are smooth functions of 1/r respecting the fall-off conditions ([5]). 
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3 The Weyl tensor 



In this section, we determine the asymptotic fall off of Weyl tensor components for asymptotically 
fiat spacetimes, as defined above. We will use Bondi coordinates since this allows us to treat even 
and odd d simultaneously. We perform our calculations using the higher dimensional Geroch-Held- 
Penrose (GHP) formalism of Ref. |18j (see appendix [A] for a review). 

3.1 Expansion of metric 

We begin with the metric written in Bondi coordinates (|6l7p . From the definition of asymptotic 
flatness we have [14jl 

A:>0 fc>0 

r d/2+k ' 2^ r d/2+k ' V zu ) 

k>0 k>0 

where in all of the summations k £ Z for even d and 2k £ Z for odd d. Equation (J7J) implies that 

= for fe < d/2 — 1 (21) 
where a/" 7 is the inverse of ujjj. In particular we have, for d > 5 

t/^g = u/ j 4 3 / 2) = ^45 = o (22) 

and 

u, IJ h?] 2) =0 (d>5), ^ /J 4 5 / 2) = i^^j (d = 5), (23) 

where all indices on 4j + ^ are ra i se( i using u/^. Some of the coefficients here have special significance. 
The Bondi mass is defined as [H] 

M( u ) = - — [ AW^du (24) 

167T J s d ~ 2 

where the integral is taken over a sphere at null infinity. In vacuum it obeys the mass decrease law 

M(u) = / h { ^h^ IJ duj. (25) 

327T J gd-2 

This demonstrates that the quantity 4j wm be non-zero when gravitational radiation is present. 
Ref. [2] showed that hjj is not constrained by the asymptotic vacuum Einstein equation: it is a 
free function in the Bondi approach, just as in 4d [HE]. hjj corresponds to Bondi's "news function". 

3.2 Null frame and connection components 

We choose a null frame (£,n,m^) for the metric given by 

3 -b ( 9 lA d d \ _ _j d 



£ = m(n\ = — — , n = mm = e I — hA— C 



£ b = e B du, n b = -(dr + gAdu), m (i) b = e^dx 7 + C'du), (26) 



3 Our notation differs slightly from that of Ref. |14j . notably in the expansion coefficients of B. 
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where e, form a vielbein for the metric hu on S d ~ 2 : hjj = euejjSij. We choose this vielbein by 
using the Gram-Schmidt algorithm starting from the basis r~^e\ where the vectors e\{x) form an 
orthonormal basis for the metric ujj on S d ~ 2 . This gives an expansion in inverse powers of r (even 
d) or y/r (odd d): 

eu = r L + ®\ + 0(r-W), ef = r" 1 U - jjfeL) + 0( r -^), (27) 



where 2e i(/ e« )( ^ = hf} and = u"e?}. 

Using the definition of the connection components given in appendix lAl and the null frame given 
in (]26p . we find that the GHP covariant connection components are 

m = 0, pij = -\e\ej d r (r 2 hu) , p = -(d-2)/r, n = -\ (a ■ dB + e~ B eud r C I ) 
k- = \e~ B ei ■ dA, p' {j = -e~ B e( i \ K \e j) ■ dC K + \e\e]n-d (r 2 h u ) , 

T '.=-^ &i .QB- e - B e lI d r C I ) . (28) 

The non-covariant coefficients are 

L\o = —d r B, Lu = ^e~ B d r A, Lu = —r[ 

Mjo = eyp r e^i, Mji = e~ ■ dC e^ K - e^n ■ de^ 

M jk = r 2 e( e/ef d {I h J]K - ef^ ■ de^. (29) 

Using the asymptotic behaviour of the metric components given in (|20p gives 

£ e 7 e J / h {1) h {3/2) \ 



„d/2 ' > r (d+l)/2 

™J 9r 2 J/2-1 T r (i-l)/2 ^j) r d-2^^v /> 



r . = 4 £^ + d+l g - c(3/2) + 0(r -w 2+ D) 

* 4 r d/2 ^ 4 r (d+i)/2 ^^V' )■> 

^-f^-^^ + o^y ,30, 

where a dot denotes a partial derivative with respect to u. Also, 

L w = ^ + + + ( r -(^)) 

2r d / 2+1 2 r ( d + 3 )/ 2 

d-2AW d-1 A( 3 / 2 ) ml M/2 ,,v r rf/2N 

8 i = r- 1 e i -V + 0(r- d / 2 ), (31) 

where V/ denotes the covariant derivative induced by w/j. Of course, terms with half-odd- integer 
powers appear only for odd d. 
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3.3 Parallely transported frame 

The null basis introduced above is convenient for calculations but it is not parallely transported 
along the geodesies. A parallely transported basis is one for which, in addition to the geodesic 

i 

equation Kj = 0, we have t[ = L\q = Mjo = 0. Any such basis will be related to ours by a boost, 
spin and null rotation (see Appendix). Let (I, n, rhi) be such a basis. I must be parallel to I, with the 
coefficient fixed by requiring that I correspond to affine parameterization of the geodesies, ensuring 
fci = I/io = 0. This gives I = e~ B £, corresponding to a boost with parameter e~ B . t[ is invariant 
under a boost and transforms covariantly under a spin. But under a null rotation with parameters 
Zi it transforms inhomogeneously [18J so Zj is determined by t[ = 0. This gives 

Zi = Ci + 0(r~W 2 -V) (32) 

i 

where the parameters c, are independent of r. Finally, Mjo transforms homogeneously under a boost 

i 

and trivially under a null rotation but inhomogeneously under a spin. Requiring Maq= determines 
the spin matrix to be 

Xij = Oij + 0(r~W 2 -V) (33) 

where 0{j is a r-independent orthogonal matrix. 

Our strategy will be determine curvature components in the basis defined previously and then 
transform our results to a parallely transported frame by first performing a boost with parameter 
e~ B = 1 + 0(r~ d / 2 ), then a null rotation with parameters Z{ and finally a spin with parameters Xij 
as given above. 

3.4 Calculation of curvature components 

In the GHP formalism the Weyl tensor components are denoted 

&ij = CoiOj, ^ijk = Coijk, V&j = CoiOi = ^jij, 

&ijkl = Cijkh ®ij = Coilj, <E> = &u = CoiOl, 

(2$| = 2* (ii) = -$ ikjk , = 2 $ m = C 01ij ), 

^ij = Cmj, ty' ijk = Cujk, *t = Clou = ^jij (34) 

and the Ricci tensor components are 

W = -RoO) = Roii 4>ij = Riji 4> = Rail 4>'i = Rlii ^ = R\\. (35) 

The Newman-Penrose equations (see Appendix) are used to determine all of these quantities except 
for those of boost weight zero (i.e. those written with the letters <3? or </>). To determine the latter 
we used the Bianchi equation (B3) from the Appendix^ 

3.5 Results: even d 

In our basis (I26p . we find that the Ricci tensor components are smooth functions of 1/r with 

u = 0(r-W 2+ V), ^ = 0(r-W 2+ V), fa = 0(r-^ 2+1 )), 

( j ) = 0{r- { - d l 2+ V), ip' i = 0(r~ d / 2 ), uj' = 0{r- d l 2 ). (36) 



4 This involves an integration with respect to r, introducing a homogeneous term decaying as 1/r into the boost 
weight zero quantities. This is not compatible with asymptotic flatness so the coefficient of this term must vanish. 



This could be shown e.g. by using the commutator (C3) of Ref. [18] . 
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The Weyl tensor components are smooth functions of 1/r with 



(d-2)(d-4) gfefog _ (d/2+2) , 



* ijk = 0(r-W 2+ V), % = <D(r-W 2+ V), ^ = 0(r-W 2+1 )), $ = 0( r -W 2+1 )), 

^5 _ (j ~ 4 ) e i e j^IJ , fl, -(d/2+lh *..„ _ (&TsJx„._ filpJ X„.\!!lJ_ +0( r -(d/2+l)) 

^ ifc = 0(r- d / 2 ), *i = 0(r- d / 2 ), 
1 p 7 p J /j {1) 

«« = -57*# + 0('- d/2 )- 

Recall that /ijj is non-zero in any spacetime containing outgoing gravitational radiation, and it is 
not determined by the asymptotic Einstein equation. 

Now we transform to a parallely transported frame as determined above. The boost and spin 
are easy to deal with since the curvature components transform covariantly with respect to these. 
Formulae for the transformation under a null rotation are given in the Appendix. Using these results, 
we see that the transformation to a parallely transported frame does not change any of these results 
(aside from acting with the rotation matrix Oij on the indices i,j,k etc). 

Finally we have to convert from our parameter r to an affine parameter along the geodesies. 
Denote the latter by A. Then (up to the freedom to multiply by a quantity independent of r) 

A = J e B dr = r + c + 0{r^ d / 2 -^) (38) 

where c is independent of r. Inverting gives 

r = \-c + (D(\- W2 - 1) ). (39) 

If we substitute this into the above expressions for the Weyl components then they become smooth 
functions of 1/A with leading order behaviour given by replacing r with A in these expressions. 
Hence the leading order term in the Weyl tensor is of order \~(- d l 2 ~ l ) and the only non-vanishing 
components at this order are D,'^, which (from (|37p ) is generically non-zero. But this is precisely the 
definition of a type N Weyl tensor with t (the tangent to the geodesies) a multiple Weyl Aligned 
Null Direction (WAND) [T7]. 

The next non- vanishing terms in the Weyl arise at order \~ d l 2 . Such terms can arise from Q^ -, 
^'ijki ^'ii ^»jW an d &fj- So, at this order, we have Q{j = ^>ijk = VPj = and hence the Weyl tensor 
is type II with multiple WAND £. It cannot be type III because ([37|) shows that $5?- is generically 
non-vanishing. However, it is not the most general possible type II Weyl tensor because (I37j) shows 
that it has vanishing $ and $4. 

After this, we have terms of order \~( d / 2+1 ) . At this order, any of the Weyl components can be 
non-zero. In particular, the above expression shows that Qij is generically non-zero, which implies 
that the Weyl tensor at this order is type G (i.e. I is not a WAND). 

In summary, for even d > 4, we have demonstrated that, in a spacetime satisfying the definition 
of asymptotic flatness at null infinity of Refs. 0[14], ^ ne Weyl tensor exhibits the peeling behaviour 
described around equation ([2]), with the type II part obeying the additional conditions $ = $4 = 0. 

When d = 4, our results for O^-, vE^ -j, and ^ are consistent with the 4d peeling property. The 
boost- weight zero terms also are consistent: in 4d, all such terms are determined by $ and which 
vanish at order \~ d / 2 = A -2 . Hence at order A -2 we have a Weyl tensor of type III instead of type 
II. More explicitly, in 4d, &ijki is determined by its trace But the first term in the expansion 
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of $>fj in (I37p comes with a coefficient of d — 4. Similar results hold for the other Weyl components 
(e.g. the above expression for has a factor c£ — 4). This is why peeling is qualitatively different 
when d = 4. 

3.6 Results: odd d 

As discussed above, for odd d there is an additional condition in the definition of asymptotic flatness 
at future null infinity, that the term of order r~( d / 2-1 / 2 ) in the expansion of hjj should be absent, 
i.e., 

4j /2) = 0. (40) 

With this condition, we find that the results (136|37p for the Ricci and Weyl components are valid also 
for odd d, with the understanding that these formulae now refer to expansions in inverse powers of 
*Jr. (Without (|40p . there would be e.g. a term of order r~( d ' 2 ~V 2 ) in the expansion of f^ -.) Following 
the same steps as for even d, converting to a parallely transported frame and affine parameterization, 
we find, just as before, that the leading components of the Weyl tensor arise at order \~W 2 ~tf and 
this term is type N as before. There are no terms at order ,\ — w 2— V 2 ) S o the next term is at order 
\~ d / 2 which is type II with <I>4 = <I> = 0, again as for even d. 

A difference between even and odd d arises at next order: for odd d there is the possibility of 
terms of order r -(<V 2 + 1 / 2 ). Yoi example, we find that 

_ l r -(d/2+i/2) faffi/*) _ i(3/2) 5 .. _ Tem ^ . V c^ K ) + 0(r-W 2+1 )), 

(41) 

where Yij is a quantity whose explicit form we will not need. The Weyl components O^- are obtained 
by taking the traceless part of this equation and the Ricci component to' by taking the trace. We 
have retained a term of order r~( d ~ 2 ^ because r - ^ -2 -* = r~( d / 2+1 / 2 ) if d = 5. 

The only other Weyl components containing terms of order r -( d / 2 + 1 / 2 ) are v[r' anc l v[r^ which 
can be obtained from 

V'm + ^^i = r-^efe^^V^j+O^-W 2 ^)), 

v^__L_^ = ^/^ ; . cW + M r -WHi/^ . (j(3/2) + o^-t^+i)) (42) 

where the Weyl and Ricci components can be disentangled by taking a trace of the first equation 
and combining with the second equation. 

We will now argue that terms of order r _ ( d / 2+1 / 2 ' can be eliminated for d > 5 by exploiting the 
Einstein equation (which we did not use for even d). We will assume that the Ricci tensor (and hence 
the energy-momentum tensor) decays faster near infinity than the rate which is given by asymptotic 
flatness alone (equation (|36p ). The rate that we require is faster by a factor 1/rH 

w = ( r -W 2 +3)) ) ^ = 0{r -(d/2+2)^ ^ = (r-W 2+ V), 

= C>(r-W 2+2 )), ^ = 0(r-W 2+1 )), to' = 0(r-W 2+ V). (43) 

For d > 5, the constraints on uj, ipi and <f>, necessarily imply the constraints on and to' . For d = 5, the constraint 
on ipi implies the constraint on f/ij. 
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Imposing these conditions implies that the first few coefficients in the expansions of the metric 
components must satisfy the same equations as in a vacuum spacetime, as determined in Ref. |14j : 



B« = 0, A® = --^-V • C« = --A-V'V- 7 ^ 1 ], C^ 1 = ~VjhW IJ , 

CI £ CLytL^ij ill 

r i( 3 /2) = i^d^w d = 5 f-hhn»h$ d = 5 

, / s (l ,A'L^(1) 1,(1) J _ r 

c^=o, hf 2) = r h Kn h J)L d - 5 (44) 

and an equation relating hj) to A^\ C^ 1 and /i^j. Note that the asymptotic Einstein equation 

implies no restriction on h^J. Recall that for d = 5, A^"/ 2 ^ determines the Bondi mass via (|24p Fl 

Using these results, we see that the term of order r~( d l 2+1 l 2 ^ in (|42p is absent and hence such 
terms do not appear in and However, terms of this order are absent from (|41h if, and only 
if, d > 5. Hence, for d > 5, such terms are absent from Q'^. Transforming to a parallely transported 
frame and affine parameterization, similar arguments to those used for the even d case establish 
the peeling result given in equation ([2]) for odd d > 5. As for even d, the type II term obeys the 
additional restrictions = $ = 0. 

Finally we must discuss the d = 5 case. For d = 5, terms of order r -(<V 2 +V 2 ) = r ~ 3 j Q no ^ ^rop 
out of f2'„-: 



0'. = _U_J__ Ll + 



•7/2) _ 
(45) 



Note that the coefficient of r~ 3 is quadratic in h^J and its time derivatives and hence generically it 
is non-zero if gravitational radiation is present. 

Now we must transform to a parallelly transported frame. As for d > 5, the boost and null 
rotation do not change our results. But note that the spin matrix (|33|) involves a term of order 
r -3/2_ fj ence w hen the spin acts on Q'- this term will combine with the leading term in O^- to 

produce a new term of order r~ 3 in Cl'-. Could this new term cancel the terms already present? 

Generically no: the new term will involve fyfj whereas some of the terms already present involve 

only first derivatives of h^J . Since h^J is a free function in the Bondi approach, these terms will not 

cancel in general. For example, one could choose hqj to be zero somewhere, with hqj non-zero. 



The last step is to convert to affine parameterization using (|39p , which does not change anything. 
We conclude that for d = 5, the Weyl tensor satisfies the peeling property ([3]) described in the 
introduction. Again the type II term obeys the additional restrictions <3?^ = <3? = 0. 



6 Note also that h^j — for d > 5 implies that one can impose the additional boundary condition h^J 2 ^ = for 
d > 5. Ref. 14] examined the vacuum Einstein equations to higher order and the results suggest that the definition of 
asymptotic flatness for odd d should be augmented with the condition = for k — 1/2, 3/2, . . . , d/2 — 2 although 

we will not assume any more than ((401) . 
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4 Bondi flux 



M{u) = — lim — / / ^i{u, r, x)du 
r-»oo Air J S 2 J_ oq 



In 4d, the rate of decrease of the Bondi energy at future null infinity is given in terms of the 
Newman-Penrose Weyl scalar ^4 as 

du (46) 

Is 2 ' 

where dui is the volume element on a unit S 2 . In d > 4 dimensions, the rate of decrease of the Bondi 
energy at future null infinity is given by (|25p [14] . We can rewrite this in terms of fi£ • (the analogue 

of ^4) as follows. Assume that the Bondi flux vanishes in the far past, i.e. hfj — > as u — > —00. 
Then from (137)) (which holds for even or odd d) we have 

/u 
n' lj (u,r,x)du (47) 
-00 

and hence 

M(u) = - lim T —— [ (f nUu,r,x)du) doo (48) 



r— >oo 



8vr 



S d-2 



where dco is the volume element on a unit S d 2 and (Yij) 2 = YijYij. In practice, the RHS is computed 
by choosing coordinates so that the asymptotic metric takes the form 

ds 2 ~ -du 2 - 2dudr + r 2 dw 2 . (49) 

One then chooses a null vector field n that approaches ±(d/du — ^d/dr) asymptotically (the sign 
does not matter here) and a set of orthonormal spacelike vectors (i = 2, ... d — 1) such that 
n ■ mty = 0. Then Q'^ = C abcd n a m b ^n c rn d ^y 



5 Discussion 



We derived our result using Bondi coordinates since this allows us to treat even and odd d together 
for much of the analysis. However, for even d it would be more elegant to derive the peeling property 
using the conformal approach. It would be nice to see this worked out. 

For odd d > 5, our result ([2]) involves only inverse half-odd-integer powers of A. Inverse integer 
powers will appear if one continues to higher orders in the expansion. It would be interesting to know 
at what order inverse integer powers first appear. If one strengthens the definition of asymptotic 
flatness as suggested in footnote [6] then it seems likely that the first such terms will appear at order 
j n a g reemen t w ith our result for d = 5. 

Ref. [19] studied asymptotically flat solutions in d > 4 dimensions that are algebraically special. 
It was found that the latter condition is incompatible with gravitational radiation (in contrast with 
the d = 4 case). We can see a similar result here: if £ is a WAND then Qij must vanish. For 

d > 4, (|37p then requires nfj = 0, which implies vanishing Bondi energy flux, i.e., no gravitational 
radiation. 
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A Higher dimensional GHP formalism 

In this appendix, we review the higher dimensional GHP formalism of |18j . Given a background 
solution, we choose a null frame (£,n,m^) such that in this frame, the metric takes the form 

9\w = 2^n^) + mj^m^. (50) 

In the GHP formalism, one breaks complete covariance by singling out two null directions {£ and 
n) at each point, but preserves covariance in the remaining directions. This is in contrast to the NP 
formalism where none of the covariance is preserved. 

At any point, the Lorentz group is generated by 

• boosts (/i a real function): 

£—>fj,£, n — > ^ _1 n, mu\ — > mu\, (51) 

• spins (Xij 6 SO(d — 2)): 

£ — > £, n — > n, — > Xijmu\, (52) 

• null rotations about £ [%i d — 2 real functions): 

£ — )■ £, n — > n + Zim^ — \z 2 £, — > — Zi£, (53) 

where A ^ and Xy is some position-dependent orthogonal matrix. 

We would like to keep the subgroup that preserves the null directions, i.e. the subgroup given 
by boosts and spatial rotations, or spins. Thus, we would like to work with objects that transform 
covariantly under this subgroup. 

A GHP scalar of boost weight b and spin s is a scalar r\i x ...i a that transforms covariantly as 

Vh...i s -> /'V....,, (54) 

under boosts and 

Vii...i a Aj ljri • • • Xj, s j s rjj lM , M j s (55) 

under spins. Evidently, the product of two GHP scalars of boost weights b\ and 62 and spins s\ and 
S2, respectively, gives a GHP scalar of boost weight b\ + 62 and spin s\ + si- 
Denote the covariant derivatives of the basis vectors as 

i 

L^v = V v £^ = V^n^, Mfj,u = V v m^ , (56) 

Not all the scalars formed from the projection of these objects into the basis are GHP scalars. Those 
that are GHP scalars are listed in table Q] [18 . 

Notice that we have used a prime operation, which interchanges the null basis vectors 

' : t^n. (57) 

For a GHP scalar t^...^ of boost weight b and spin s, we define its GHP covariant derivatives to be3 

k 



\>Ti 1 i 2 ...i s — £ ■ dTi 1 i 2 ...i a ^10^112. + Mi r oTi 1 ...i r -iki r +i---i s > (58) 

r=l 
S k 

r=l 

^i^hh—js = m (i) ' dTj 1 j 2 ...j s — bLuTj 1 j 2 ^j s + Mj r jTj 1 ___j T _ 1 kj r+1 ^j s . (60) 



7 Symbols b and 6, pronounced "thorn" and "eth", respectively are old Germanic letters that have been retained in 
the Icelandic alphabet. 
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Table 1: GHP scalars constructed from covariant derivatives of the basis vectors. 



Spin coefficient 


GHP notation 


Boost weight b 


Spin s 


Interpretation 


La 


O ; a 


1 


2 


expansion, shear and. twist of £ 


La 




1 


o 


pxnansion of / 






2 


1 


non-geodesity of t 


Ln 


n 





1 


transport of I along n 


Nij 


p'ij 


-1 


2 


expansion, shear and twist of n 


Nu 


p' = P'a 


-1 





expansion of n 


N a 


< 


-2 


1 


non-geodesity of n 


N i0 


< 





1 


transport of n along I 



GHP versions of the " Newman-Penrose" and Bianchi equations are given in Ref. [18J. We will 
need the NP equations: 

\>Pij-bjKi = -pikPkj ~ KiTj - TiKj - fiy - j=2U)6ij, (NP1) 

pri-p'Ki = Pij (-Tj + r' j )-% + - s ^tPi, (NP2) 

2 %\Pi\k\ = 2r iP[jk] + 2K iP[jk] - ^ijk ~ d=2^[j 5 k]i, (NP3) 

Vpij-Sjn = -nTj - KiKj - p ik p' kj - Qij - + <t>Sij) + ( / fc i fc )ff- 2) ^i- (NP4) 

Another four equations can be obtained by taking the prime ' of these four (i.e. by exchanging the 
vectors £ and n). We also need the following components of the Bianchi identity: 

- \>$ ijkl + 26 [fe % j = -2* ' m K {j] - 2& [klij K {l] 

+^tjP[kl] ~ 2®[k\iPj\l] + 2 ${k\jPi\l] + 2®ij[k\mPm\l] 

- 2 ^[i\kl T '\j] ~ 2 ^[k\ij T '\l] ~ 2 ^i[k\p'j\i] + 2 ^j[kP'i\i], (B3) 

where the tilde notation indicates components of the Riemann tensor analogously defined to those of 
the components of the Weyl tensor given in (|34p . The relation of these components to the components 
of the Weyl and Ricci tensors is given by the definition of the Weyl tensor as the trace-free part of 
the Riemann tensor 

W r T 2 / r t 1 / 

;Oij, ^ijk = ^ijk ~ ~, — ~W\jOk\h ^i = ^i + -, — -m 







&ijki 


= ®ijkl 




i 




d-2 




= % ~ 



d _ 2 ->->> -v» -v». d _ 2 T w^J" ' d-2* : 

2 / , „ , „ x 2 



d 



i[k$l]j - < Pj[k^l]i) + 7^ _ l^d - 2) ^ + ( t }mm ^i[k^l]j-> 



{(pdij + 0ijj + J- 7T7^ ^T^i) * = * + 



(d-l)(d-2) d-1 (d-l)(d-2)' 

*^ = %fc-^2^%> *i = *i + ^2$- (61) 
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Null rotations 



Under a null rotation about i of the form given by equation (|53p the Weyl tensor components 



transform as: 

^ + (63) 

^ijk !-»■ + (64) 

$ i-> $ + 2zi^i + ZiQijZj, (65) 

*ij ^ $y +^*i + ^i^+%Oi fc , (66) 

$yfcj !->■ - 2^[ fe *j]y - 2^*^ - 2zjZ[ k Qi]i + %ZiZ[ktti]j, (67) 

¥j ^ tf^ - + 3$Azj- - $gz,- - 2Z ij V j - Z jk y jki - ZjZ ik Q jk , (68) 

+2Z u z [k n j]b (69) 

^y ^ fi y ~ 2z (J% + 2z fc*(i|*W + 2Z m®k\j) + - Az k z (i <5>f )k + z k zi$ mj 

+2z( i Zj) k ^ k + 2ziZ{i\ k ^ k i\^ + Z ik ZjiVt k i. (70) 
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